Modification of surface energy in nuclear multifragmentation. 
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Within the statistical multifragmentation model we study modifications of the surface and sym- 
metry energy of primary fragments in the freeze-out volume. The ALADIN experimental data on 
multifragmentation obtained in reactions induced by high-energy projectiles with different neutron 
richness are analyzed. We have extracted the isospin dependence of the surface energy coefficient at 
different degrees of fragmentation. We conclude that the surface energy of hot fragments produced 
in multifragmentation reactions differs from the values extracted for isolated nuclei at low excita- 
tion. At high fragment multiplicity, it becomes nearly independent of the neutron content of the 
fragments. 

PACS numbers: 25.70.Pq , 25.70.Mn , 21.65.+f 



I. INTRODUCTION 

A break-up of nuclei into many fragments (multifrag- 
mentation) has been observed in nearly all types of nu- 
clear reactions when a large amount of energy is de- 
posited in nuclei. It includes reactions induced by pro- 
tons, pions, antiprotons, and by heavy ions of both, 
relativistic energies (peripheral collisions) and 'Fermi'- 
energies (central collisions) . According to the present un- 
derstanding, multifragmentation is a relatively fast pro- 
cess, with a characteristic time around 100 fm/c, where, 
nevertheless, a high degree of equilibration is reached. 
The process is mainly associated with abundant produc- 
tion of intermediate mass fragments (IMFs, with charges 
Z « 3-20). However, at the onset of multifragmentation, 
also heavy residues are produced which have previously 
only been associated with compound-nucleus processes. 
At very high excitation energies, the IMF production 
gives way to the total vaporization of nuclei into nucleons 
and very light clusters. 

The previous ALADIN experiments have provided ex- 
tensive information about multifragmentation of projec- 
tiles in peripheral nucleus-nucleus collisions at high en- 
ergy 0. They have demonstrated a 'rise-and-fall' of mul- 
tifragmentation with excitation energy, and they have 
shown that the temperature remains nearly constant, 
around T ~5 MeV, during this process 0- There are 
large fluctuations of the fragment multiplicity and of the 
size of the largest fragment in the transition region from a 
compound-like decay to full multifragmentation of spec- 
tators jli. It was found that the statistical models, which 
assume a thermal equilibration among hot fragments in a 
freeze-out volume at subnuclear densities, are fully con- 
sistent with the data |1I1,H,|3|- 

We believe that multifragmentation studies pursue two 
main purposes. The first one is connected with the gen- 
eral understanding and better description of this reaction 



which represents as much as 10-15% of the total cross sec- 
tion in high energy hadron- nucleus collisions, and nearly 
twice more in nucleus-nucleus collisions. Secondly, the 
multifragmentation reaction can be considered as an ex- 
perimental tool to study the properties of hot fragments 
and the phase diagram of nuclear matter at densities 
p « 0.1 — 0.3/90 and temperatures around T w 3-8 MeV 
which are expected to be reached in the freeze-out vol- 
ume (po 0.15 fm^^ is the normal nuclear density). 
Multifragmentation opens a unique possibility to investi- 
gate this part of the phase diagram and to determine the 
"in- medium" modifications of nuclei there. This second 
point is very important for many astrophysical applica- 
tions, in particular, for processes duringthe supernova II 
explosions and neutron star formation |3, 0) ■ 

Some of the properties of hot fragments have been ad- 
dressed in our previous work. For example, the symme- 
try energy was extracted in ref. and it was demon- 
strated that it decreases considerably with excitation en- 
ergy. In this paper we show that also the surface energy 
of hot fragments can be investigated in these reactions, 
and that the modified surface energy can be extracted by 
comparing theory with experiment. 



II. STATISTICAL APPROACH FOR THE 
DESCRIPTION OF MULTIFRAGMENTATION 

Statistical models are used in situations when an equili- 
brated source can be defined in the nuclear reaction. The 
most famous example of such a source is the 'compound 
nucleus' introduced by Niels Bohr in 1936. The standard 
compound nucleus picture is valid only at low excitation 
energies when sequential evaporation of light particles 
and fission are the dominating decay channels. However, 
this concept cannot be directly applied at high excita- 
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tion energies, E* > 2-3 MeV/nucleon, when the nucleus 
rapidly disintegrates into many fragments. As was shown 
in many experiments (see e.g. 0,0,0,0,0,0]), an 
equilibrated source can be formed in this case too, and 
statistical models are generally very successful in describ- 
ing the fragment production. 

As a basis for our study below we take the Statisti- 
cal Multifragmentation Model (SMM), for a review see 
ref. 0. The model assumes nuclear statistical equi- 
librium at a low-density freeze-out stage. It considers 
all break-up channels composed of nucleons and excited 
fragments taking into account the conservation of baryon 
number, electric charge and total energy. Light clusters 
with mass number A < 4 are treated as elementary parti- 
cles with only translational degrees of freedom (" nuclear 
gas" ) . Fragments with A > 4 are treated as heated liquid 
drops. In this way one may study the nuclear liquid-gas 
coexistence in the freeze-out volume. The Coulomb inter- 
action of fragments is described within the Wigner-Seitz 
approximation. Different decay channels / are generated 
by Monte Carlo sampling according to their statistical 
weights, oc expSf, where 5/ is the entropy of channel 
/. After the break-up, the Coulomb acceleration and 
the secondary de-excitation of primary hot fragments are 
taken into account. The SMM has been already success- 
fully apphed for the analysis of ALADIN data [l3,0|- 
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III. INFLUENCE OF THE SYMMETRY 
ENERGY ON MULTIFRAGMENTATION 



In the SMM, the symmetry energy of hot fragments 
with mass number A and charge Z is parametrized as 
^a'z ~ ~ 2Z)^/yl, where 7 is a phenomenological 
parameter. In the case of cold nuclei 7 w 25 MeV is 
adopted in order to describe their binding energies. For 
hot fragments this parameter can be modified and, there- 
fore, should be extracted from experimental data. The 
corresponding set of isospin-related data may be provided 
by the multifragmentation reactions, e.g., via the isoscal- 
ing phenomenon |0|. At present, there are evidences 
for a significant reduction of the symmetry energy in hot 
nuclei [lli^. 

The influence of the symmetry energy on the character- 
istics of produced fragments was already partly analyzed 
in [23. It was shown that it has a small effect on av- 
erage quantities such as the temperature and fragment 
masses, but it influences mainly the isotope distributions 
of the fragments. This is illustrated in Fig. 1 for the mean 
charge distributions and neutron-to-proton (N/Z) ratios 
of fragments. These characteristics change only slightly 
even with rather large variations of 7. In this respect, 
we believe that variations of 7 will not influence the re- 
sults of previous analyses of experimental data related to 
the fragment char ge p artitions and average isotope yields 
0000000111. 

By contrast, the symmetry energy has a large effect 
on isotope distributions. As evident from Fig. 2, they 
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FIG. 1: Influence of the symmetry energy coefficient 7 on 
yields of liot fragments in the freeze-out volume (bottom pan- 
els) and on their N/Z ratios (top panels), in SMM calculations 
for Au sources at excitation energies of 3 and 5 MeV/nucleon 
(left and right panels, respectively). 



become wider at small 7 |19l |22| . The subsequent deex- 
citation of the produced hot fragments changes the dis- 
tributions in two ways: their widths usually become nar- 
rower and their centroids are shifted toward symmetry 
[T^ I21I I23I I. This effect has to be taken into account 
when comparing with experimental data: it is essential, 
e.g., for the isoscaling analysis pdL Il9l l2(]||. 

In the following, we study possible modifications of 
the fragment surface energy employing different isospin 
compositions of the sources. Within the SMM this issue 
can be separated from modifications of the symmetry en- 
ergy because the latter does not significantly change the 
mean charge yields. In phenomenological mass formulae, 
isospin dependent terms for the surface energy are some- 
times introduced in order to obtain a better description of 
the ground-state masses of nuclei [3|2l. However, some 
modern mass formulae can provide even better descrip- 
tions of the nuclear masses retaining only the bulk isospin 
term but with a special treatment of the shell effects z6] . 
We believe that the best strategy for the present purpose 
of investigating new 'in-medium' properties of nuclear 
fragments is to consider a minimum number of param- 
eters which have a clear physical meaning. We, there- 
fore, use the standard SMM liquid-drop parametrization 
which contains no surface isospin term. Instead we allow 
for the surface coefficient to be a free parameter which 
is determined by comparing the calculations with exper- 
imental data obtained for multifragmentation of nuclei 
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FIG. 2: Mass distributions of hot fragments with atomic num- 
bers Z — 6, 12, and 26, produced in the multifragmentation of 
Au sources at excitation energies of 3, 5 and 10 MeV/nucleon, 
for different symmetry energy coefficients 7. 



FIG. 3: Influence of the surface energy coefficient Bq on 
charge yields of hot fragments in the freeze-out volume, as 
obtained from SMM calculations for Au sources at excitation 
energies of 3, 5, and 10 MeV/nucleon. 



with different isospin content. 



IV. INFLUENCE OF THE SURFACE ENERGY 
ON MULTIFRAGMENTATION 

The surface free energy of hot fragments in the SMM 
is parametrized as F^'J = B{T)A^/^, where B{T) = 
Bo[(r2 - T^)/{T^ + r2)]5/4. Here Bq «18 MeV is the 
surface coefficient for isolated cold nuclei, as adopted 
in the Weizsaecker mass formula, and Tc ~ 18 MeV is 
the critical temperature for the nuclear liquid-gas phase 
transition in infinite matter. One should distinguish Tc 
from the phase transition temperature in finite hot nu- 
clei, which is essentially lower, around 5-6 MeV i5|. Tc 
should be considered as a model parameter characteriz- 
ing the temperature dependence of the surface tension 
in finite nuclei. At low temperatures this T-dependence 
leads to a correct surface contribution to the level densi- 
ties of nuclei, i.e., it describes correctly the heat capaci- 
ties of isolated nuclei. The surface parameters may also 
be modified in the low-density nuclear medium, i.e., in 
environments consisting of nucleons and hot fragments. 
The possible changes of in multifragmentation reac- 
tions were analyzed in refs. |^ |23|, |28|, and from the 
analysis of experimental data it was found that Tc ~ 17- 



20 MeV in reactions with Au nuclei. The Tc effect on 
fragment yields is rather small and, in any case, it is 
included in the surface energy. In the present work we 
concentrate on the influence of the Bq coefficient on the 
multifragmentation pattern. 

The surface energy is rather important, since produc- 
tion of new fragments means increasing the surface con- 
tribution to the total energy of the system. Therefore, 
even small variations of the surface energy lead to big 
changes of fragment mass and charge distributions. This 
is demonstrated in Fig. 3 for different excitation ener- 
gies of Au sources with a freeze-out density p = po/3. 
Decreasing the surface energy favors the disintegration 
already at low excitation energies while a larger Bq sup- 
presses the multifragmentation channels. 

To characterize the charge distributions we use the 
Z~'^ fit of the fragment yields. In order to avoid con- 
tributions of other processes leading, e.g., to fission- like 
large fragments {Z > 20) and to coalescence-like small 
clusters {Z<2), we limit the fragment charges by the 
range Z = 5 — 15 in this fit. The extracted r parame- 
ters for three sources with different isospin, ^^®U, -'^^^Au 
and ^^^Xe, are shown in Fig. 4. In these calculations, 
the surface coefficient was fixed at the standard value 
of Bq = 18 MeV. One can see that the isospin effect is 
rather essential, as was already discussed in [2lll27| . 



4 




4 6 8 10 
E* (MeV/n) 

FIG. 4: SMM calculations of r parameters for hot fragments 
from •^'^^U, '^■'^Au and '^'^'^Xe sources as a function of their 
excitation energies. The power-law fitting wEis performed for 
fragments with Z = 5 — 15. 



We choose the Au source to demonstrate possible ef- 
fects of the modified surface energy on multifragmenta- 
tion. Figure 5 shows the influence of the coefficient Bq on 
the caloric curve, on the mass number Amax of the largest 
fragment, and on r. We should point out that thermo- 
dynamical characteristics of the system in the freeze-out 
volume can also be influenced by the modifications of the 
surface energy of fragments. When these hot fragments 
leave the freeze-out volume and decay, their normal prop- 
erties are restored. In order to emphasize the difference 
between isolated nuclei and nuclei in the medium, we in- 
troduce two different temperatures. One is an effective 
temperature, Teg, which is found from the energy balance 
in the freeze-out volume by assuming that the properties 
of fragments are the same as those of isolated nuclei. 
This effective temperature reflects internal excitations of 
fragments respective to their ground states and, thus, 
it can be compared to the temperature reached in the 
compound nucleus at the same excitation. Another tem- 
perature is the freeze-out temperature T which includes 
in- medium modifications of the fragment properties, i.e., 
different Bq. The evolution of both temperatures with 
E* is presented in Fig. 5 in the top two panels. One can 
see that variations in Bq lead to noticeable effects. With 
a small Bq, correlated with a smaller effective tempera- 
ture, the system disintegrates into lighter fragments, and 
^max becomes smaller too. At low excitations the effec- 
tive and freeze-out temperatures exhibit a similar behav- 
ior. However, at very high excitation energies when the 
system disintegrates only into light IMFs, the freeze-out 
temperature becomes higher at smaller Bq. This behav- 
ior has a simple explanation: the smaller surface term 
leads to a smaller level density parameter for internal 
excitation of fragments, which requires higher freeze-out 
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FIG. 5: SMM calculations of characteristics of hot frag- 
ments from Au sources for different surface energy coefficients 
Bo = 14, 18, 22 MeV as a function of the excitation energy 
E*. Top two panels: effective temperature Tgg, and freeze- 
out temperature (caloric curves); middle panel: reduced mass 
number of the maximum fragment Amax/^o; bottom panel: 
T parameter. 



temperatures in order to achieve the energy balance. The 
T parameter exhibits a general trend in the evolution of 
the charge distributions, it depends strongly on the ex- 
citation energy for all considered values of the surface 
coefficient (Fig. 5, lower panel). 

The change in the disintegration process can be con- 
veniently characterized by Tmin, i-e- the minimum value 
which r assumes as a function of the excitation energy. 
The region near Tmin corresponds to the transition from 
channels with one big residual fragment to multifragmen- 
tation into several small fragments. In the nuclear matter 
this transition is usually associated with the coexistence 
region of the liquid-gas phase transition. Both, Tmin and 
the corresponding effective temperature Tmin, vary as a 
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function of the surface energy, as demonstrated in Fig. 6. 
As one can see, the three considered sources with differ- 
ent isospin content exhibit similar trends but the r,„in for 
the neutron-rich sources is higher. 
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FIG. 6: Minimum value Tmin of the r parameter for hot frag- 
ments and the corresponding effective temperature Tmin, ex- 
tracted from SMM calculations for sources with different neu- 
tron richness ^^*U, ^^'^Au, and ^'^^Xe, as a function of the 
surface energy coefficient Bo- 



V. EXPERIMENTAL DATA ON THE 
T-PARAMETRIZATION OF THE CHARGE 
DISTRIBUTIONS 

In order to extract information about the surface en- 
ergy in muhifragmentation we compare the calculations 
with experimental data. The ALADIN collaboration has 
performed exclusive measurements of muhifragmentation 
events for projectile-like sources with different isospins 
and has made these data available. Measurements have 
been made with ^'^*U, ^^^Au, and -"^^^Xe projectiles at the 
laboratory energy of 600 MeV per nucleon, interacting 
with different targets ranging from Be to U. The power 
law parameters r, obtained from fits to the charge distri- 
butions in the range 5 < Z < 15, are shown in Fig. 7 as 
a function of the bound charge Abound (the total charge 
accumulated in fragments with Z > 2) divided by the 



projectile charge Zq. Small 2'bound values correspond to 
high excitation energies of the sources which disintegrate 
predominantly into very light clusters ("fall" of multi- 
fragmentation). Large values of .Z^bound correspond to 
low excitation energies, at which the decay changes its 
character from evaporation/fission like processes to mul- 
tifragmentation ("rise" of multifragmentation). 

The chosen target has no influence on the extracted t 
parameter, as demonstrated in the bottom part of Fig. 7 
for the case of the Au and U projectiles. This is consis- 
tent with the target invariance generally observed for the 
charge correlations of projectile fragmentation , !]• These 
invariance properties extend also to projectiles with dif- 
ferent Zq if the data are scaled accordingly. For example, 
the IMF multiplicities divided by Zq (the reduced multi- 
plicities) are practically identical if plotted as a function 
of the reduced bound charge Zbound/^o- This has been 
shown in for the case of U, Au and Xe projectiles with 
the same data sets as used here. The same holds for the 
T parameters which, however, exhibit some variations at 
small and large .Abound (Fig- 7, top). Differences at small 
■Abound may occur for several reasons. The uncertainty 
of determining Abound for highly fragmented partitions, 
arising from the finite acceptance of the spectrometer 
for He and Li fragments may be slightly system de- 
pendent which, in principle, could be taken into account 
by using an experimental filter. Restricting 2'bound to 
small values will, furthermore, generate constraints for 
the accepted partitions which affect the resulting Z dis- 
tributions 3] . The fitting within a fixed Z interval while 
scaling Zbound could thus possibly lead to systematic de- 
viations in T for different projectiles which, however, are 
small as evident from the data. 

At large Abound, the observed differences are slightly 
larger and more systematic, and no experimental con- 
straints of this type exist. The detection efficiency of 
the ALADIN spectrometer is, furthermore, very high for 
Z > 5, essentially An in the projectile frame, and the 
experimental trigger was adapted for the fragmentation 
channels T]. In the following, we will, therefore, concen- 
trate on the analysis of the observed differences in the be- 
havior of T for sources with different isospin at the "rise" 
of multifragmentation (.^boundZ-^o ^ 0.5). Because of the 
large acceptance of the detector, the theoretical predic- 
tions will be directly compared to the data without the 
use of an experimental filter. 



VI. THEORETICAL CALCULATIONS WITH 
ENSEMBLES OF EQUILIBRATED SOURCES 

In the most general consideration the fragmentation 
process can be subdivided into several stages as fol- 
lows. (1) A dynamical stage leading to the formation 
of equilibrated nuclear systems, (2) disassembly of the 
system into individual primary fragments, and (3) de- 
excitation of hot primary fragments. Even though the 
first stage may be described by a dynamical model, such 
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FIG. 7; Power law parameters r obtained from fitting tlie 
fragment cliarge distributions obtained in ALADIN experi- 
ments in the range 5 < Z < 15 as a function of the re- 
duced bound charge Z\,ound/Zo for the fragmentation of ^'^^Xe, 
^'^''Au and ^^**U projectiles of 600 MeV/nucleon. Top panel: 
three systems with different projectiles as indicated; bottom 
panel: results for collisions of ^^^Au with Be, C, Al, Cu, In, 
Au targets and of ^''^U with Cu, In, Au, U targets. 



as the intranuclear cascade model, it is more practical 
to determine an ensemble of equilibrated sources, pro- 
duced after the nonequilibrium stage, by direct com- 
parison with experimental data 0|. In the following, 
we consider the ensembles of equilibrated sources ex- 
tracted in refs. 0,0] for the Au projectiles, which pro- 
vide a very good description of the ALADIN data. We 
use also the same version of the SMM code described 
fully in 0. As was found in ref. the average 

masses of equilibrated sources As may be parametrized 
as As/Ao = 1 - ax{E*/As) - a2{E*/Asf, where E* is 
the excitation energy of the sources in MeV, and 
is the projectile mass. Here is ai=0.001 MeV^^, and 
02 «0. 009-0. 015 MeV~^ is slightly depending on the pro- 
jectile energy. It should be noted that the adopted corre- 
lation of the mass number with the excitation energy is 
consistent with dynamical simulations |2 
with results of other statistical models 



sol, as well as 



[Q. Besides 

providing a very good description of fragment produc- 
tion at different .Abound, the calculations with this kind 
of ensemble reproduce also correctly the behavior of the 
caloric curve [ig. 

Our present analysis is moderately sensitive to the rela- 
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FIG. 8: Ensemble of hot thermal sources represented in a 
scatter plot of reduced mass number As/Ao versus excitation 
energy E* /As for the fragmentation of ^^^Au projectiles, as 
used in the SMM calculations. The intensity of the individual 
sources is proportional to the area of the squares. 



tion between As and E* because we extract t in different 
Zhonad/Zo bins, which are correlated with the excitation 
energy (per nucleon) of the source. Possible variations of 
the relative weights of the sources in the ensembles have 
practically no effect on the r observable, as was confirmed 
by performing calculations with different weights. The 
typical relative yields of sources used for the case of ^^''Au 
projectiles are shown in Fig. 8 as a function of their mass 
and excitation energy. The parameter a2=0.012 was used 
in this case. The same relation between masses and exci- 
tation energies is also taken for the other projectiles. It is 
assumed that the N /Z ratio of all sources in the ensemble 
is the same as in the projectiles. The dynamical calcu- 
lations show that this is a quite reasonable assumption 
for sources with moderate excitation energies. A very 
small decrease of the N/Z ratio is predicted only for the 
sources with very high excitation energy pjl Il9| | . This 
should not qualitatively influence the present analysis in 
the region of the multifragmentation rise. 

The SMM calculations with these ensembles for differ- 
ent projectiles predict a scaling of the fragment multiplic- 
ities, as well as other observables, in agreement with the 
experimental data a. This is illustrated in Fig. 9, where 
we plot the reduced IMF multiplicities versus ^bound/^o- 
In order to demonstrate the sensitivity of this correlation 
to the ensemble parameters and to the surface energy, we 
compare calculations obtained with different parameter 
sets. For -'^^^Au we show results for the ensemble pa- 
rameters 02=0.011 (full squares) and a2=0.012 (empty 
squares), obtained with the surface term i3o=18 MeV. 
For 238u we compare the cases Bq=\1 MeV (empty trian- 
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FIG. 9: Mean IMF multiplicities versus ^bound/^o, the bound 
charge divided by the projectile charge, predicted by the SMM 
ensemble calculations for three diflterent projectiles. The 
meaning of the symbols is explained in the text. 



FIG. 10: SMM ensemble calculations of r parameters for cold 
fragments versus Zhound/Zo for ^^*U, ^^'^Au, and ^'^^Xe pro- 
jectiles. The surface energy coefficient is -Bo=18 MeV for all 
sources. 



gles) and Bq—IS MeV (full triangles) for ensembles gen- 
erated with 02=0.011 and, for the ensembles representing 
the fragmentation of ^^^Xe, generated with 02=0.012, we 
show the results for Bo=19 MeV (full circles) and Bo—20 
MeV (empty circles). In this way we cover the range of 
parameters representative for the following analysis of the 
data. One can clearly see a very good scaling behaviour, 
similar to the experimental one demonstrated in Fig. 9 
of ref. f|. A good description of the IMF multiplicities 
allows us to adjust the ensemble parameters and start 
the T analysis. This analysis allows us to extract more 
detailed information on the multifragmentation process, 
since the r fit in the range Z=5-15 represents a new sen- 
sitive observable. It is important to mention that the 
maximum IMF production is reached at E* around 5-6 
MeV/nucleon, when the sizes of the sources are relatively 
large (Fig. 8). 

The results of the SMM calculations for ensembles of 
the sources with a2=0.012 produced with Xe, Au, and 
U projectiles are shown in Fig. 10. In the case of Au, 
the obtained t parameters are close to the experimen- 
tal values. This remains also valid in the region of the 
multifragmentation "fall" , i.e. at high excitation energies 
and small .Abound, even though the sources become small 
and have very large size variations (Fig. 8). However, 
in the region of the multifragmentation "rise" the the- 
ory gives larger r for the neutron-rich U projectile than 
for the neutron-poor Xe projectile. This is opposite to 
the experimental observation. Below we investigate pos- 
sible modifications of the SMM which can explain this 
discrepancy. 

First, we have tested whether one can improve the 



agreement with the experiment by simply using mass 
formulae with an explicit isospin dependence of the sur- 
face energy. We have implemented the Myers-Swiatecki 
mass formula in the SMM, instead of the stan- 
dard liquid-drop description of hot fragments. In this 
case the surface energy coefficient is written as Bq = 
18.56(1 - 1.79(1 - 2Z/Af) MeV. The other parame- 
ters were also modified according to ref. The re- 
sults of the SMM calculations with the Myers-Swiatecki 
parametrization are shown in Fig. 11 for the ensembles 
with 02=0.012, for all three projectiles. The resulting 
values of r for different projectiles at large Abound are 
less spread than the results with constant Bq — 18 MeV 
(Fig. 10) but they do not fit the experimental trend. 
Moreover, the general agreement with the experiment 
becomes even worse, since the minimum of r is consid- 
erably shifted to larger values of .Abound- Also Tmin is 
less than 1.5, which is too low in comparison with the 
experimental value of about 1.7. Calculations with other 
ensembles did not improve the agreement with the exper- 
iment. This suggests that a more detailed investigation 
of the liquid-drop parameters of fragments in multifrag- 
mentation reactions is needed. 

Returning to the standard SMM, we have performed 
calculations for the same projectile ensembles but with 
other surface coefficients Bq. We have selected four in- 
tervals Zbound/^o=0.8-0.9, 0.7-0.8, 0.6-0.7, and 0.5-0.6, 
characterizing the "rise" of multifragmentation with ex- 
citation energy, for a detailed comparison with the ex- 
periment. As one can expect from Fig. 5, in order to fit 
the higher r observed for neutron-poor Xe sources at low 
excitation energies, it is necessary to increase the sur- 
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FIG. 11: SMM ensemble calculations of r parameters for cold 
fragments versus Zbound/^^o for '^'^^U, ^^'^Au, and ^'^^Xe pro- 
jectiles. The Myers-Swiatecki mass formula was used for the 
hot primary fragments. 



face coefficient Bq for the corresponding hot fragments. 
As known from previous analyses (see, e.g., |27,]), the 
secondary decay processes should not significantly influ- 
ence T in this region. For the secondary deexcitation 
calculations we have used the standard evaporation and 
Fermi break-up models implemented in the SM M J23| . 
The same deexcitation code was used in refs. [UlUfor 
describing the ALADIN data. 

The theoretical evolution of r with Bq in different 
•Abound bins is shown in Fig. 12. In order to achieve 
agreement with the experimental data, Bq should change 
with excitation energy when going from compound-like 
processes to full multifragmentation. In the bin of the 
largest Zhound/Zo = 0.85, this requires Bq « 20.0, 18.1, 
and 17.1 MeV for the Xe, Au, and U ensemble sources, 
respectively. At Zbound/^o=0.65 and 0.55, however, in 
the region of minimum r and full multifragmentation, 
the values of the surface tension are around Bq «18-19 
MeV for all three sources without essential variation. 

We have also analyzed other ensembles and combina- 
tions of ensembles for different projectiles, which can 
reproduce the IMF scaling shown in Fig. 9 and in ref. 
[l'|. In addition, we require that these ensembles repro- 
duce the experimental t in the whole range of .Abound by 
varying Bq. We have found that such ensembles should 
have parameters in between a2=0.011 and a2=0.013. For 
control, we have also tested the extreme cases as, e.g., 
a2=0.011 for the ensembles of U sources and a2=0.013 
for the Xe sources. Even for these ensembles, we have 
obtained a similar evolution of the Bq isospin dependence 
as a function of the excitation energy: while Bq depends 
strongly on N/Z at large 2^bound, it becomes nearly in- 
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FIG. 12: Comparison of the r parameters obtained from the 
SMM ensemble calculations as a function of Bo (dashed lines 
and open symbols) with the experimental results (solid hori- 
zontal lines and full symbols) for U, Au, and Xe projectiles, 
in different Abound /^o bins. 



dependent of isospin in the region of full multifragmen- 
tation. 



Additional calculations were performed with another 
secondary decay code developed in ref. This code 

starts with the modified surface tension for hot fragments 
in the freeze-out volume, and then changes the fragment 
masses during the evaporation process in such a way that 
the normal surface tension is restored for cold fragments. 
It was found that the final results may deviate from the 
standard calculations by not more than a few percent, 
and that all trends regarding Bq as a function of the 
neutron richness of fragments remain the same. 
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VII. DISCUSSION OF THE RESULTS 

As was postulated earlier fSect. lTVj) . the surface term is 
a function of two parameters, Bq and the critical temper- 
ature Tc- Therefore, a dependence of on the neutron 
content could be an alternative explanation of the isospin 
dependence of the surface tension in different domains of 
excitation energy. This possibility was investigated by 
performing an analysis with a possible Tc evolution, as 
suggested in refs. [2^ |2^. It was found that, in order 
to increase r to the experimental value observed for Xe 
sources at large Zbound, 7c would have to be increased 
to Tc « 24 MeV. In the same Abound range, in order to 
reproduce the value for U sources, Tc would have to be 
decreased to Tc « 14 MeV. A Tc variation of this mag- 
nitude in a narrow N/Z range (from 1.39 to 1.59) is not 
supported by theoretical calculations |3ll |. which predict 
a variation of Tc within 2-3% only. In addition, large 
Tc differences for the projectile sources will result in a 
significant violation of the scaling behaviour of fragmen- 
tation at high excitation energies, in contradiction to the 
experiment. 

We have mentioned that for small sources, at 
Zbound/^o <0.5, the experimental filter should be ap- 
plied for a precise determination of Abound hi simulated 
events. This would be important for our method since at 
high excitation energy, as one can see from Figs. 4 and 
10, the r extracted for the sources of different isospin 
nearly coincide, while the r values after the ensemble cal- 
culations are different, when plotted versus the reduced 
bound charge. This is a consequence of dealing with the 
small sources and large fluctuations at the multifragmen- 
tation 'fair. Also, in the region of small .Abound, after 
SMM calculations, many observables are less sensitive to 
differences in the ensembles of sources: As evident from 
Fig. 9, the slope of Mj^jf/Zq versus ^bound/^o at small 
Zbound is rather universal for all ensembles and values of 
Bq, and it coincides with the data 0|. Without indepen- 
dent constraints on the sources in this region, it is not 
very reliable to extract fragment properties by only using 
the T parametrization of the charge yields. Probably, in- 
volving isotope distributions of produced fragments and 
measuring neutrons j34l 13^ will help to determine the 
ensembles. We do not touch this problem here but just 
note that extraction of Bq as a function of the neutron 
richness at small Abound will depend essentially on the 
excitation energies of the sources. As one can see from 
Fig. 5 (bottom panel), if the average excitation energy 
is around 7-8 MeV/nucleon, the sensitivity of r to the 
surface coefficient is lost at Bq « 18 — 22 MeV. 

Besides their different N/Z ratios, the three projectiles 
differ also in their mass. To first order, this is considered 
by performing the analysis as a function of scaled vari- 
ables but a residual mass effect may possibly exist. The 
calculations presented in ref. show, however, that 
for relatively large sources Ag ~ 100 a difference in the 
N/Z ratio has a larger effect on r than a difference in 
the mass. Moreover, calculations for the fragmentation 
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FIG. 13: The extracted Bo for ensemble sources with different 
N/Z ratios at the onset (Zbound/.Zo=0.85 bin), and at the 
region of full multifragmentation (Zbound/^o=0.65 and 0.55 
bins). The two dotted lines for the lowest bin indicate upper 
and lower limits given by our method. The Bo obtained from 
Cameron (C) and Myers-Swiatecki (MS) mass formulae for 
cold nuclei are shown by solid lines for illustration. 



of nuclei with different N/Z and the same mass predict 
differences in t of the same magnitude as obtained here 
for the considered Xe, Au and U nuclei (Fig. 10). This 
justifies our conclusion that the deviation of the stan- 
dard SMM calculations from the experimental data at 
low excitation energies is due to the isospin dependence 
of the surface coefficient. This conclusion can be tested 
with fragmentation data in experiments with isotopically 
different isobars. 

It is instructive to compare the surface coefficients ex- 
tracted above with the ones used usually for cold iso- 
lated nuclei. In Fig. 13 we show the behavior of this 
coefficient versus the N/Z ratio as extracted from the 
Cameron and Myers-Swiatecki [2^ mass formulae 
by retaining terms proportional to A"^'^. The consid- 
ered N/Z range corresponds to the neutron-richness of 
the sources under investigation. Both mass formulae 
predict a decrease of Bq with increasing N/Z which is 
necessary in order to reach agreement with experimen- 
tal mass tables. In multifragmentation, as was shown in 
refs. 2Qi ..22] , the neutron content of hot fragments 
is proportional to the neutron richness of the sources, 
since most neutrons are accumulated in the fragments. 
The values of Bq extracted from the experiment exhibit 
the same decreasing trend with N/Z at low excitation 
energy and they remain within the Bq range predicted 
by the mass formulae. Thus, at the rise of multifrag- 
mentation, when a heavy residue and few small frag- 
ments are formed, Bq resembles the values used in phe- 
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nomenological parameterizations for isolated nuclei. On 
the other hand, at full multifragmentation, when the sys- 
tem breaks up into intermediate mass fragments, the sur- 
face energy coefficient becomes nearly independent of the 
neutron richness. This is particularly evident for the bins 
^bound/^o=0.65, where the decrease in Bq with N/Z is 
very small, and Zbound/^o=0.55, at which the Bq consis- 
tent with T falls between the two dotted lines (Fig. 13). 
These upper and lower limits correspond to the 2% sta- 
tistical uncertainty in the numerical determination of r, 
typical for both, the SMM simulations and the experi- 
mental data. 



These observations indicate that nuclear fragments are 
brought to a new physical environment. Their proper- 
ties change in the hot medium consisting of nucleons 
and other fragments at subnuclear freeze-out densities, 
around O.lpo- Also, modifications of the syrnmetry en- 
ergy have been observed in these reactions [ll|. Here, we 
do not suggest a comprehensive theoretical explanation 
of these effects, but rather point at the new physical con- 
ditions reached in multifragmentation reactions in com- 
parison with isolated nuclei. The origin of these effects 
could be a residual nuclear interaction between clusters, 
in particular, the exchange of neutrons between nuclear 
fragments. The role of the Coulomb forces is also impor- 
tant. The Coulomb interaction between clusters leads to 
a modification of their average charges |2^. Other con- 
sequences of the Coulomb interaction, like a shift of the 
proton distributions with respect to the neutron distribu- 
tions in hot fragments js^, may additionally contribute 
to this effect. 



Finally we point to the uncertainty connected with the 
lacking knowledge of the N/Z ratio of the produced equi- 
librated sources after the initial dynamical stage. The 
present dy namical codes, like the intranuclear cascade 
(INC) [13, and the relativistic quantum molecular dy- 
namics (RQMD) 0,1^, predict a very small, less than 
5%, change in comparison with the neutron richness of 
the projectile nucleus, since interactions of neutrons and 
protons are very similar at relativistic energies. This 
change is expected to be small for all projectiles, there- 
fore, the difference in neutron content should be pre- 
served for the ensembles of projectile sources. Similar 
conclusions can also be drawn on the basis of existing 
experimental data: As was established in recent AL- 
ADIN experiments, fragments coming from neutron-rich 
projectiles remain essentially neutron-rich as c omp ared 
with fragments from neutron-poor projectiles [S^. At 
the same time, the number of produced spectator neu- 
trons is correlated directly with the neutron richness of 
the projectiles [s^. These results indicate that the ini- 
tial difference in N/Z ratios for different projectiles is 
not washed out for the ensemble of sources during the 
nonequilibrium stage. 



VIII. CONCLUSIONS 

In this paper we have demonstrated that the surface 
properties of hot fragments can be studied in multifrag- 
mentation reactions. Within the Statistical Multifrag- 
mentation Model we have shown how modifications of 
the fragment surface energy influence the fragment pro- 
duction and thermodynamical characteristics of the sys- 
tem at the freeze-out volume. By comparing with the 
ALADIN data we have found that in the regime of the 
" rise" of multifragmentation the observed isospin depen- 
dence of the power-law parameter r of fragment charge 
distributions cannot be described by the SMM if the stan- 
dard surface-term coefficient is used. The data can be 
reproduced if a moderate dependence on the N/Z ratio 
of the produced fragments is introduced in the surface 
term. At low excitation energies this contribution re- 
sembles that obtained with advanced mass formulae for 
ground-state nuclei. At higher excitation energies, this 
N/Z dependence gradually disappears. 

The observed modification of properties of hot nuclei in 
comparison with cold (or slightly excited) isolated nuclei 
are possibly caused by the hot environment in which the 
nuclei are imbedded, and where they can interact with 
nucleons and other fragments. The resulting in-medium 
modifications can go beyond the disappearance of shell 
effects and affect also the main liquid-drop parameters. 
Our conclusions are substantiated by a very good agree- 
ment with experimental data demonstrated by the SMM. 
In the future, fully microscopic calculations of the effects 
of the interacting medium would be very desirable, how- 
ever, many problems still remain in constructing such a 
many-body theory. 

More generally, we conclude that multifragmentation 
reactions offer a new possibility for investigating nuclei 
under conditions essentially different from those accessi- 
ble in nuclear structure studies at low energies. Besides 
the obvious application of the results for better describing 
nuclear reactions at intermediate and high energies, this 
has far-reaching consequences for astrophysics. In many 
astrophysical situations, hot matter at subnuclear den- 
sities is clusterized into nuclei as, e.g., in neutron star 
crusts and supernovae. The properties of these nuclei 
are crucial for both microscopic (weak interactions) and 
macroscopic (collapse dynamics and shock propagation) 
processes in stellar matter. New and important informa- 
tion in this direction may be expected from the analysis 
of future experimental data on multifragmentation. 
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